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•^ . Abstract 

o . 

^b , We introduce a very simple, exactly solvable T-'T-symmetric non-Hermit- 

f~^ ' ian model with real spectrum, and derive a closed formula for the metric 

\^ \ operator which relates the problem to a Hermitian one. 



1 Introduction 



9^: 

rl^ ' In a way motivated by the needs of nuclear physics, Scholtz, Geyer and Hahne PP 

j^ , established a general mathematical framework for the consistent formulation of 

quantum mechanics where a set of observables are represented by bounded non- 
Hermitian operators Ai^ . . . , An with real spectra in a Hilbert space Ti. In 
essence, they conjectured that in the similar situations one may find a bounded 

^^ , positive Hermitian operator 0, called metric, which fulfils 

C3 ; AlQ^QAk foraU fc e {1, . . . , TV} , (1) 

where A*j, denotes the adjoint operator of Ak in TL. 

Several years later, the notion of the metric operator 8 re-emerged as a par- 
ticularly useful mathematical tool in the context of the so-called VT -symmetric 
quantum mechanics [3 El- In this framework people usually paid attention to 
the systems with a single observable, viz, with a Hamiltonian Ai = H ^ H* 
which possesses real spectrum and for which the Schrodinger equation is invari- 
ant under a simultaneous change of spatial reflection V and time reversal T. 

In the current literature a lot of effort has been devoted to the study of 
the particular models of H . For their more detailed reviews and discussion 
the reader is referred to the proceedings of the International Workshops on 
Pseudo-Hermitian Hamiltonians in Quantum Physics 0] |S1 . One finds that 
the construction of a non-trivial operator O =/= I, however difficult, is a key to the 



correct probabilistic interpretation of all the PT-symmetric quantum systems 
0111111 En]. Indeed, it defines "the physical" inner product {■,-)e := (-.O-) 
which makes the Hamiltonian H "Hermitian" or, in the language of 0, quasi- 
Hermitian. For this reason, there have been many attempts to calculate the 
metric operator Q for the various PT-symmetric models of interest m El El 
^^ El El El El El ■ Because of the complexity of the problem, however, it is 
not surprising that most of the available formulae for are just approximative, 
usually expressed as leading terms of perturbation series ^71 . 

The authors of PP discussed why our knowledge of the new inner product 
was necessary for the evaluation of the physical predictions. They emphasized 
that the theory endowed with it is a genuine quantum theory satisfying all 
the necessary postulates. In a fairly recent continuation of this discussion PU] 
it has been underlined that in the infinite-dimensional Hilbert spaces H. the 
requirement of the boundedness of the metric operator Q plays a key role and 
that it deserves an extremely careful analysis in the applications where a naive 
approach may lead to wrong results. In some sense, our present paper may be 
read as a direct continuation of the rigorous mathematical discussion in ^E^ . 

In particular we are going to illustrate here that our understanding of (QJ for 
unbounded operators H as the identity on functions from the operator domain 
oi H {cf lO below) requires that 9 maps the operator domain of H into the 
operator domain of the adjoint H* . In such a setting we imagined that the best 
way of finding a support for such an argument can be sought in some exactly 
solvable "PT-symmetric model. We decided to develop a new one - such that 
its metric can be obtained in a closed formula and in a rigorous manner. 

The model we deal with in the present paper is one-dimensional, defined in 
the Hilbert space 

H:=L^{iO,d)) 

where d is a given positive number. In this Hilbert space we consider the 
Hamiltonian H^ which acts as the Laplacian, i.e., 

H^^ := -^" , 

and satisfies the following Robin boundary conditions: 

^'{0) + iaip{0)^0 and ip'{d) + ia^^d) = . (2) 

Here i/' is a function from the Sobolev space W^'^{{0, d)) and a is a real constant. 
That is, the operator domain D{Ha) consists of functions with integrable (gen- 
eralized) derivatives up to the second order and satisfying Q at the boundary 
points. Because of the nature of the boundary conditions, H^ is not Hermitian 
unless a = 0, but it is "PT-symmctric with the operators V and T being defined 
by (Vtp)ix) :— 'ip{d — x) and Tip := -0, respectively. 

It seems that our Hamiltonian Ha represents the simplest PT-symmetric 
model whatsoever. The fact that the support of the non-Hermitian perturba- 
tion is of measure zero invokes the "PT-symmetric models [21 1221 E] involving 
complex point interactions. But our model is even simpler, since it does not 



require any matching of solutions known explicitly off the points where the 
(5-interaction is supported. 

Indeed, the non-Hermiticity of Ha enters through the boundary conditions 
only, while the Hamiltonian models a free quantum particle in the interval (0, d). 
Consequently, the spectral problem for Ha can be solved explicitly in terms of 
sines and cosines (c/ Section |31 for more details). Furthermore, an explicit form 
for the eigenfunctions enables us to obtain a remarkably simple expression for 
the metric operator: 

Theorem 1. Let Q{a) he the linear operator defined in TL by 

e(a) := / + 0;? (C ■) + 60 + la^i + a^Oa , (3) 

where I denotes the identity operator in Ti, (■,•) is the inner product on Ti, 
antilinear in the first factor and linear in the second one, 

<Pq{x):^ J-ey.^{iax) (4) 
and the operators Oq, Qi and 02 acts in Ti as 

(eoV)(a:):=-i(J^)(d), (5) 

(eiV')(x) := 2 {J^){x) - ^ {JrP){d) - i (jV)(d) , (6) 

{Q2i^){x) := -{J'^){x) + 5 (^V^)(rf) , (7) 

with 

(J^)(x):=rV. (8) 

Jo 

Then 0(a) is bounded, symmetric, non-negative and satisfies 

y^eD{Ha), H*aeia)ij^eia)Hai'. (9) 

Furthermore, 0(a) is positive if the condition 

ad/ir ^ Z\{0} (10) 

holds true. 

Note that the metric 0(a) tends to the identity operator / as a ^ 0, which is 
expected due to the fact that Hq is nothing else than the (self-adjoint) Neumann 
Laplacian in TC. The condition H1U|) ensures that all the eigenvalues of Ha are 
simple. For simplicity, we do not consider the degenerate cases in the present 
paper. 

This paper is organized as follows. In the following Section |21wc introduce 
the Hamiltonian Ha by means of its associated quadratic form; this provides 
an elegant way of showing that the operator is closed. The spectral problem 



for Ha is considered in Sectional in particular, we show that the spectrum is 
real and discrete, and write down the explicit eigenfunctions and eigenvalues. 
Section ^contains the main idea of the present paper; namely, we observe that 
the eigenfunctions of Ha are expressed in terms of Dirichlet and Neumann com- 
plete orthonormal families in the interval (0, d) and use a special normalization 
to simplify the eigenfunctions of the adjoint 7J*. These enable us, in Section |31 
to evaluate certain infinite series defining the metric operator and prove Theo- 
rem ^ We conclude the paper by Section El where we add several remarks and 
discuss a possible extension of our model. 

2 The Hamiltonian 

Let us first introduce the operator H^, in a proper way. We start with the 
associated sesquilinear form h^ defined in the Hilbert space Ti, by the domain 
D{ha) ■= W^''^{{0,d)) and by the prescription: 



ha{(j), ij) := (</)', ij') + ia (t>{d) ^{d) - ia (/.(O) ^(0) . (11) 

Here (•,•) denotes the standard inner product on Ti; the corresponding norm 
will be denoted by || ■ ||. 

Note that the boundary terms in Hll|) are well defined because the domain 
of the quadratic form is embedded in the space of uniformly continuous func- 
tions on (0, d) due to the Sobolev embedding theorem j21]. It is also known 
that the Sobolev space W^''^{{0,d)) is dense in Ti; hence ha is densely defined. 
Moreover, the real part of ha, denoted by ^ha, is a densely defined, symmet- 
ric, positive, closed sesquilinear form (since it corresponds to the self-adjoint 
Neumann Laplacian in Ti). Of course, ha itself is not symmetric unless a = 0, 
however, it can be shown that it is sectorial and closed. To see it, we use [251 
Thm. VI. 1.33] and prove that the imaginary part of ha, denoted by '^ha, is a 
small perturbation of "Siha in the following sense: 

Lemma 1. Qha is relatively bounded with respect to ^ha, with 

\{^haM\ <e-^a^U\\^ + e{maM 
for all ip £ W^''^{{0, d)) and any positive constant e. 

Proof. Writing |^(d)P - |i/'(0)|2 = /„'' (1^1^)' = 23?(V',V'')' ^nd applying the 
Schwarz and Cauchy inequalities to the last term, we obtain the desired result. 

D 

In view of the above properties and the first representation theorem |25l 
Thm. VI. 2.1], there exists a unique m-sectorial operator Ha in H such that 
ha{4>,ii^) = {(j),Haii^) for all cf) G D{ha) and i/" € D{Ha) C D{ha)- The operator 
domain D{Ha) consists of those functions i/" G D{ha) for which there exists 



r] G 7i such that ha{4>,ip) — {4',v) holds for every </> G D{ha)- Furthermore, 
using the ideas of '5S!, Ex. VI. 2. 16], it is possible to verify that indeed 



ip e D{Ha) ^{^& M^2'^((0,d)) I V satisfies ©} 



2,2,.. ,^^, , ...... ™. (12) 



The above procedure also implies that the adjoint operator 7J* is simply ob- 
tained by the replacement a t^ —a. 
Summing up the results, we obtain: 

Proposition 1. The operator Ha defined by ]1}^) is m-sectorial in Ti and sat- 
isfies 

3 The spectrum 

An important property of an operator being ?7i-sectorial is that it is closed. 
Then, in particular, the spectrum is well defined by means of the resolvent 
operator. We claim that our H^ is an operator with compact resolvent. This 
can be seen by noticing that the Neumann Laplacian Hq (associated with Wia) 
is an operator with compact resolvent and by using the perturbation result 
of |25[ Thm. VI. 3. 4] together with Lemma Q] Consequently, we know that the 
spectrum of iJ^, denoted by a{Ha), is purely discrete, i.e., it consists entirely 
of isolated eigenvalues with finite (algebraic) multiplicities. 

The eigenvalue problem Haip — k'^ip, with fc G C, can be solved explicitly in 
terms of sines and cosines. In particular, the boundary conditions lead to the 
following implicit equation for the eigenvalues: 

{k^ - a^) sm{kd) ^ . (13) 

That is, 

a{Ha) = {«'} U {fcfl^'li , where k, := jir/d . (14) 

Hereafter we shall use the index j e N to count the eigenvalues as in (|14|l , with 
the convention that the eigenvalue for j = is given by a^ . 

While the spectrum of Ha is real, it exhibits important differences with 
respect to the spectra of self-adjoint one-dimensional differential operators. For 
instance, the spectrum of Ha may not be simple and even the lowest eigenvalue 
may be degenerate for particular choices of a. Notice also that Ha coincides with 
the spectrum of the Neumann Laplacian Ho up to the lowest (zero) eigenvalue 
which is shifted to a^. 

In this paper we restrict to the non-degenerate case, i.e., we make the hy- 
pothesis (|1(J|) . Then the eigenfunctions of Ha corresponding to H14|) with the 
convention mentioned there are given by 

\ Aq exp (—iax) if j = , 

^j"(^) •= ] ^a (r.r.^(i,.^\ „ ^»_ e;„,^.^^^ if ,■ -> 1 (1^) 



1" (cos(fcjx) — i-^ sin(fcjx)] if j > 1 



where each Aj is an arbitrary non-zero complex number. In view of Proposi- 
tion^ the spectrum of the adjoint if* coincides with 1)141) and the corresponding 
eigenfunctions are given by 

{Bg exp (iax) if j = , 

Bf (cos(fcjx) + if- sm{kjx)] if j > 1 , ^^^^ 

where each Bj is again an arbitrary non-zero complex number. 

We collect the obtained spectral results into the following proposition: 

Proposition 2. The spectrum of Ha is real and consists of discrete eigenvalues 
specified in ^14^ . If the condition ^^ holds, then all the eigenvalues have 
multiplicity one and the corresponding eigenfunctions are given by ^^. 

4 Special normalization 

It follows directly by combining the eigenvalue problems for Ha and its adjoint 
that 0" and ipk ^^^ orthogonal to each other provided j ^ k and the non- 
degeneracy condition (|10(l holds. The stronger result 

Vj-fceN, {^'-,i,<^)^5,k (17) 

will hold provided we normalize the eigenfunctions in a special way. Namely, 
l(T7|) follows by choosing the coefficients A" and i?" according to the equations 

l = ^„"Wi^^^^|^, (18) 

ik^-a^)d 

1 = ^j Bf \^2 fo-^ J > 1 • (19) 

(If a = 0, the fraction in the first equation should be understood as the expres- 
sion obtained after taking the limit a -^ 0.) These equations can clearly be 
satisfied as soon as ((1(1)1 holds. 

We still have a freedom in specifying A" and B". For further purposes, 
however, we choose the coefficients B" in a very simple form by the requirements 

Bo := ^/T/d and Bj := y^ for j > 1 , (20) 

while we leave more complicated formula, determined by the equations H18|l 
and H19I) . for the coefficients A". Then 0q coincides with Q) and we have the 
decomposition 

<f>f{x)^xf{x)+t^xfix) for j->l, (21) 

Kj 

where {xj^lj^O' respectively {xf}^i^ denotes the set of normalized eigenfunc- 
tions of the Neumann, respectively Dirichlet, Laplacian in Ti: 

Xf(x):=|^ ^, ^ '^, ^Z^' xf(x):-y^sin(fc,x). 



In addition to H21I) . we also have the uniform convergence (f>Q -^ Xo as a ^ 0. 
We point out the result we shall need later: 

Proposition 3. // the condition \1U(I holds true, then the eigenf unctions tp" 
of Ha and the eigenf unctions (j)" of _ff * can be chosen in such a way that they 
satisfy the biorthonormality relations \11^ and the latter satisfy ^21\) . 

The decomposition H21II plays a crucial role in the subsequent section, mainly 
due to the fact that {xf}?Lo a-nd {xf }?°.i are well known to form complete 
orthonormal families. In particular, we have the expansions 

oo oo 

^ = E ^f (Xf » = E 4 (Xf ' ^) (22) 

and the Parseval equalities 

oo oo 

uf = T.\(4,n' = T.\(xfM' (23) 

for every %jj £71. 

5 The metric 

With an abuse of notation, we initially define 

oo 

e(a):=^0,"(0,",-) (24) 

and show that this operator can be cast into the form |j2Jl with Q-Q). In fact, 
using (|21|l and 1)22(1 . it is readily seen that Q holds with 



eo:=-X^Mx6\-) (25) 



and 



°° yD ( N \ _ N ( D \ oo D ( D \ 

ei:^J2 ' k ' 92 :^ J2 I2 ' ■ (26) 

Recalling the definition (jHJ of the bounded integral operator J in ?i, it is evident 
that the rank-one operator (|25|l can be expressed in terms of J as in (0. It 
remains to verify that H2()|) can be expressed as in © and {7)). 

First of all, we notice that the operator 1)24(1 is well defined in the sense 
that Oi and 82 as defined in l(2()() are bounded linear operators in H. This can 
be seen by using ((23(1 and the Schwarz inequality. Actually, the series in ((26(1 
are uniformly convergent, and Q2 can be written as an integral Hilbert-Schmidt 
operator, but we will not use these facts. Our way to sum up the infinite series 
is based on the following lemma: 



Lemma 2. 



y -^ — = — — umjormly for all x G [0, d\ 



i=i 



kj d 



Proof. The series is uniformly convergent due to Abel's uniform convergence 
test. Let I denote the identity function on (0, d), i.e. l(x) := x. Using the 
expansion (|22|) and integrating by parts, we get 

oo oo oo 

^ = E ^? (^?' = E ^f ii-xf/kj)', = - E xf xf{d) d/k, , 

where the last equality follows from the fact that all Xj^ with j > I are orthog- 
onal to the constant function Xo ■ This concludes the proof. D 

Since Jtjj is an indefinite integral of 4' a-nd (J'0)(O) = 0, an integration by 
parts yields for every ip £ H: 

(xf , ^) = k, (xf , J^) + xf{d) {J^M , 

(xf , V^) = ~k, (xf , J^) = -k] (xf , jV) - % xf (d) (J'V')(d) • 

Incorporating these identities into (|26|1 and using l|22(l together with Lemma El 
we obtain the formulae © and Q for (|2SIl- 

Now we are in a position to prove Theorem ^ 

Proof of Theorem^ The boundedness of ^ is clear; in particular, crude esti- 
mates yield 

||e(a)V'|| < (3 -f 4ad + 2a^<f) U\\ 

for every ip £ TC. 

Integrating by parts, it is also easy to check that the identity 

(^, e(a)^) = !(</)?, V)l' + 11^ + ^aJ^pf ~ \ix^, ^ -f *aJ^)|2 (27) 

holds for every -0 G 7i, where the right hand side is real- valued and non-negative 
due to H23() . This proves that 0(a) is symmetric and non-negative. 

Let us show that 8(a) is positive provided l|10|l holds. If the right hand 
side of 1)27(1 were equal to zero with a non-zero tp G TC, then the first Parseval 
equality in 1)23(1 would imply that the function ip + iaJip must be constant, being 
orthogonal to all functions orthogonal to 1. Consequently, ip is proportional 
to ipQ and an explicit calculation yields 



m^n 



sin(ad) 



ad 
which is clearly positive for all a satisfying p()(l . 



Finally, let us comment on the identity ®. Let tp G D{Ha). We first note 
that it straightforward to check that ld{a)ip belongs to D{H^), so that the left 
hand side of Q makes sense. We also have 

-(e(a)i/;)" == -V/' - 2ia V' + aV + a^0o (C . V^) = -e(a)^" . 

Here the first equality follows at once, while the second one is not trivial, but 
it can be verified by using a number of integrations by parts. 

This concludes the proof of Theorem ^ D 

6 Concluding remarks 

6.1 Alternative proofs of the reality of the spectrum 

Recall that PT-symmetry itself is not sufficient to guarantee the reality of the 
spectrum of a non-Hermitian operator (see, e.g., |26[l27p . Moreover, the exist- 
ing proofs of the reality ,2gl|2niEniE| are based on rather advanced techniques. 
Therefore we find it interesting that the reality of the eigenvalues of our Hamil- 
tonian H^ can be deduced directly from the structure of the operator, without 
solving the eigenvalue problem explicitly. 

To see it, we rewrite the eigenvalue problem Haip = ^^"0 using the unitary 
transform ip i—> (jyQip := (j) into the boundary value problem 

j -(j)" + 2ia(t)' + a'^(f>^k'^<j) in (0, d) , 

I 0' = at 0, d . ^ ' 

Now we multiply the first equation in H28|l by (p" and integrate over (0, d). We 
also multiply the complex conjugation of the first equation in H28|l by 0" and 
integrate over (0, d). Then we subtract the results and use various integrations 
by parts together with the Neumann boundary conditions to get the identity 

3(fc')ll0'f = 0. 

Consequently, either the eigenvalue fc^ is real or the corresponding eigenfunc- 
tion 4> is constant. It remains to realize that also the latter implies the former 
in view of (jSHJ- 

Finally, let us mention that Ha can be reconsidered as a self-adjoint operator 
in a Krein space [29' . Then the reality of the spectrum of Ha for \a\ < n/d 
follows from [291 Corol. 3.3]. An alternative proof of the reality of the spectrum 
of Ha for small a also follows from the perturbation result of |30] . 

6.2 Biorthonormal basis 

It is easily seen that the operator 8(a) defined by (|24|l formally satisfies lO, with 
the inverse given by 9(a)"^ = Z^i^o ^f ('^'j^ ') ' provided {ip'j}JLo and {0"}°^o 



fulfil in addition to (|17|1 the following biortlionormal-basis-type relation: 

oo 

VV'GH, ij^Y.^J^'i^'j^^y (29) 

By "formally" we mean that one has to justify an interchange of summation and 
differentiation. We did not pursue this direction in the present paper. Instead, 
we summed up the infinite series (I24|l using the special normalization H2U|I lead- 
ing to (|21|l . and checked that the resulting operator indeed satisfies ^ in the 
sense of ©. 

Nevertheless, let us show that the expansion H29|l holds: 

Proposition 4. // the condition ^10\) holds true, then the eigenf unctions V'? 
of Ha and the eigenf unctions (p" of i?* can be chosen in such a way that iHyp 
is satisfied. 

Proof. Assume the special normalization of Section 0] Let us first verify that 
{V'iljLo is a basis oiTi, i.e., 

oo 

V^IjEH, V^ = Es^^"' (30) 

j=o 

where {c^}°^q is a unique sequence of complex numbers. Note that the equality 
in ()30(l should be understood as a limit in the norm topology of Ti.; in particular, 
H3(J|I implies the weak convergence 

m 

\fq>,^PeH, (0,^)= lim Uy^cf^f). (31) 

n=l 

Substituting ■)/' = and = 0^, fc G N, into (|31|l . the biorthonorniality re- 
lations p7|) yield that H3U|) with ?/; = implies that all c^- = 0. At the same 
time, 

and since the right hand side behaves as 0(j^'^) as j -^ oo, we have 

oo 

Consequently, {V'j}?=o i^ ^ basis of H due to ^', Thm. V.2.20]. Finally, substi- 
tuting (j) = (j)'^, k € N, into H31() . the biorthonorniality relations H17() yield that 
c^ = (0^,1/,) for all j eN. D 

The same argument also implies the following expansion: 

oo 
j=0 



10 



6.3 A more general model 

For simplicity, we required that a was real in the present paper. A more general 
model is given by the following more general PT-symmetric boundary condi- 
tions: 

V''(0) + (/3 + ia)V'(0) ==0 and - V^'C^) + (/? - Ja) V^d) = , (32) 

where a and /3 are real constants. A straightforward modification of the ap- 
proach of Section|21(c/ also the first paragraph of Section|2Il yields: 

Proposition 5. The operator Ha.p defined in Ti by 

V' e D{Ha^p) = {?/'£ iy2'2((o,d)) I i/, satisfies \EB] , 

is an m-sectorial operator with compact resolvent and satisfies iJ* ^ = H^a,i5- 

The eigenvalue problem Ha,p ip = k'^'^P, with fc G C, can again be solved in 
terms of sines and cosines, and one gets the following implicit equation for the 
eigenvalues: 

[k^ - {a^ + 0^)] sm{kd) -213k cos{kd) = . 

The main difference with respect to the case (3 — studied in the present paper 
is that Ha,p can possess non-real complex conjugate eigenvalues for P ^ 0. 

Acknowledgment 

One of the authors (D.K.) would like to thank the foundation Cesky literdrm 
fond for a financial support which enabled him to participate in the 3rd In- 
ternational Workshop on Pseudo-Hermitian Hamiltonians in Quantum Physics 
(Istanbul, 2005). The work has partially been supported by the Czech Academy 
of Sciences and its Grant Agency within the projects IRP AV0Z10480505 and 
A100480501 

References 

[1] Scholtz FG, Geyer HB, Hahne FJW. Quasi-Hermitian operators in quantum 
mechanics and the variational principle. Ann Phys 1992;213:74-101. 

[2] Bender CM, Boettcher S. Real spectra in non-Hermitian Hamiltonians having 
PT symmetry. J Math Phys 1998;80(24):5243-5246. 

[3] Bender CM, Boettcher S, Meisinger PN. 'PT-symmetric quantum mechanics. 
J Math Phys 1999;40(5):2201-2229. 

[4] Proceedings edited by M Znojil. 1st International Workshop on Pseudo-Hermitian 
Hamiltonians in Quantum Physics. In: Czech. J. Phys.. vol. 54; 2004. p. 1-156. 

[5] Proceedings edited by M Znojil. 2nd International Workshop on Pseudo- 
Hermitian Hamiltonians in Quantum Physics. In: Czech. J. Phys.. vol. 54; 2004. 
p. 1005-1148. 



11 



[6] Proceedings edited by M Znojil. 3rd International Workshop on Pseudo-Hermitian 
Hamiltonians in Quantum Physics. In: Czech. J. Phys.. vol. 55; 2005. p. 1045- 
1192. 

[7] Bender CM, Brody DC, Jones HF. Complex extension of quantum mechanics. 
Phys Rev Lett 2002;89:270401. 

[8] Mostafazadeh A. Pseudo-Hermiticity versus PT symmetry: The necessary con- 
dition for the reality of the spectrum of a non-Hermitian Hamiltonian. J Math 
Phys 2002;43:205-214. 

[9] Mostafazadeh A. Pseudo-Hermiticity versus PT symmetry: II. A complete char- 
acterization of non-Hermitian Hamiltonians with a real spectrum. J Math Phys 
2002;43:2814-2816. 

[10] Mostafazadeh A. Pseudo-Hermiticity versus PT symmetry: III. Equivalence of 
pseudo-Hermiticity and the presence of antilinear symmetries. J Math Phys 
2002;43:3944-3951. 

[11] Mostafazadeh A. Hilbert space structures on the solution space of Klein-Gordon- 
type evolution equations. Class Quantum Grav 2003;20:155-171. 

[12] Bender CM, Meisinger PN, Wang Q. Calculation of the hidden symmetry operator 
in ■PT-symmetric quantum mechanics. J Phys A 2003;36:1973-1983. 

[13] Bender CM, Brod J, Refig A, Renter ME. The C operator in PT-symmetric 
quantum theories. J Phys A 2004;37:10139-10165. 

[14] Bender CM, Brody DC, Jones HF. Scalar quantum field theory with a complex 
cubic interaction. Phys Rev Lett 2004;93:251601. 

[15] Bender CM, Brody DC, Jones HF. Extension of "PT-symmetric quantum mechan- 
ics to quantum field theory with cubic interaction. Phys Rev D 2004;70:025001. 

[16] Bender CM, Brody DC, Jones HF. Erratum: Extension of PT-symmetric quan- 
tum mechanics to quantum field theory with cubic interaction. Phys Rev D 
2005;71:049901(E). 

[17] Mostafazadeh A, Batal A. Physical aspects of pseudo-Hermitian and PT- 
symmetric quantum mechanics. J Phys A 2004;37:11645-11679. 

[18] Bender CM, Cavero-Pelaez I, Milton KA, Shajesh KV. PT-symmetric quantum 
electrodynamics. Phys Lett B 2005;613:97-104. 

[19] Mostafazadeh A. PT-symmetric cubic anharmonic oscillator as a physical model. 
J Phys A 2005;38:6557-6569. 

[20] Kretschmer R, Szymanowski L. Quasi-Hermiticity in infinite-dimensional Hilbert 
spaces. Phys Lett A 2004;325:112-117. 

[21] Albeverio S, Fei SM, Kurasov P. Point interactions: PT-Hermiticity and reality 
of the spectrum. Lett Math Phys 2002;59:227-242. 

[22] Fei SM. Exactly solvable many-body systems and pseudo-Hermitian point inter- 
actions. Czech J Phys 2004;54:43-49. 

[23] Znojil M, Jakubsky V. Solvability and PT-symmetry in a double- well model with 
point interactions. J Phys A 2005;38:5041-5056. 

[24] Adams RA. Sobolev Spaces. New York: Academic Press; 1975. 

[25] Kato T. Perturbation Theory for Linear Operators. Berlin: Springer- Verlag; 
1966. 



12 



[26] Znojil M. PT-symmetric square well. Phys Lett A 2001;285:7-10. 

[27] Znojil M, Levai G. Spontaneous breakdown of "PT-symmetry in the solvable 
square-well model. Mod Phys Lett A 2001;16:2273-2280. 

[28] Dorey P, Dunning C, Tateo R. Spectral equivalences, Bethe ansatz equations, and 
reality properties in "PT-symmetric quantum mechanics. J Phys A 2001;34:5679- 
5704. 

[29] Langer H, Tretter Ch. A Krein space approach to PT-symmetry. Czech J Phys 
2004;54:1113-1120. 

[30] Caliceti E, GrafB S, Sjostrand J. Spectra of PT-symmetric operators and pertur- 
bation theory. J Phys A 2005;38:185-193. 

[31] Shin KG. On the reality of the eigenvalues for a class of "PT-symmetric oscillators. 
Commun Math Phys 2002;229:543-564. 



13 



